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Abstract 

We use real replicas to investigate stability of thermodynamic homogeneity of the free energy of 
the Sherrington-Kirkpatrick (SK) model of spin glasses. Within the replica trick with the replica 
symmetric ansatz we show that the averaged free energy at low temperatures is not thermody- 
namically homogeneous. The demand of minimization of the inhomogeneity of thermodynamic 
potentials leads in a natural way to the hierarchical solution of the Parisi type. Conditions for the 
global thermodynamic homogeneity are derived and evaluated for the SK and p-spin infinite range 
models. 



I. INTRODUCTION 



Thermodynamic potentials have to be state variables. It is so if they are ther- 
modynamically homogeneous. Thermodynamic homogeneity means that the densities 
of thermodynamic potentials depend on the extensive variables only via their densi- 
ties. This property is usually expressed as the Euler condition vF{T, V,N,..., Xi, . . .) = 
F(T, uV, uN, . . . , uXi, . . .), where v is an arbitrary positive number and Xi exhaust all ex- 
tensive variables. It is straightforward to obtain from the Euler condition / = F/V = 
f(T, 1, N/V, . . . , Xi/V, . . .). This linear homogeneity condition leads to the existence of the 
thermodynamic limit independently of boundary conditions and the shape of the volume. 

One of the authors recently argued that simple solutions (replica symmetric) of the SK 
model do not obey thermodynamic homogeneity l|. Using the TAP approach a hierarchical 
solution equivalent to the Parisi replica-symmetry breaking (RSB) was derived from the 
demand of homogeneity, or better from minimizing the inhomogeneity of the averaged free 
energy. 

In this paper we show that the same demand of thermodynamic homogeneity can be 
implemented into the replica trick for the SK and other mean-field models as scale invariance 
of the replica index. 



II. REPLICA TRICK WITH REAL REPLICAS 



The replica trick is used in systems with quenched randomness to overcome averag- 
ing of logarithm. The averaged free energy is defined in the replica trick as f3F av = 
— lim ((Z n ) — 1) /n. To investigate stability of thermodynamic homogeneity we replicate 

n— >0 

v times the original spin space. The averaged free energy of such a replicated system is 
/3F„ = - (lnZ v ) av jv = - lim ({[Z u } n ) av - 1) jnv. If the limit n -> leads to a thermody- 

n—*0 

namically homogeneous solution the resulting free energy F v should be independent of the 
parameter v. Hence thermodynamic homogeneity is equivalent to scale invariance of the 
replica representation, n — > nv. 

We introduce two types of replica indices in the replica trick applied to ay times repli- 
cated system. Real replicas we denote a, b — 1, . . . , v while a, (5 = 1, . . . , n we use for the 
mathematical ones. After averaging over the random spin-spin interactions in the SK model 
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we obtain at the saddle-point 

( n v n v ~\ 
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We denoted L w = Ela 
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+ nu/2. Parameters (<3 afe ) Q/3 = (Sa^p) are deter- 



mined from saddle point equations |2|. 

To proceed further with thermal averaging we have to determine symmetry of the matrix 
(Q ab ) a/3 - We use the replica symmetric ansatz and distinguish only diagonal and off-diagonal 
elements. The super-diagonal elements are fixed, {Q aa ) aa = 1- Since the elements in the 
matrix of mathematical replicas are u x u matrices, we can choose one parameter for the off- 
diagonal elements in the real-replica space even if the indices for the mathematical replicas 
are equal. We denote (Q a ^ 6 ) aa — Qi- In the spirit of the replica-symmetric ansatz we choose 
for the off-diagonal matrix elements in the space of mathematical replicas (Q ab ) a — %■ 
For convenience we denote x = Qi ~ <?o> Q — %■ It is now a straightforward task to perform 
averaging over the spin configurations in Eq. using the above parameters. We arrive at 

h(q, x; v) = -j (i - <?) 2 + - i)x(2g + x) + fx- 

1 /*oo POO 

/ Vrjln VX{2cosh\j3(h + ri y /q + X y /x)]Y . (2) 

J —oo J — OO 



We used an abbreviation for the Gaussian differential T>(j) = d0 e - * Z 2 / \^2~k. 

The obtained free energy density depends formally on the scaling parameter v. It is clear 
from the dependence of f\ on v in Eq. (J2J) that it actually depends on v if and only if x > 0. 
We obtain easily from the saddle-point equations for q and x that this happens just below 
the de Almeida-Thouless instability line. Hence, the replica-symmetric solution becomes 
thermodynamically inhomogeneous (free energy does not obey the Euler condition) in the 
spin-glass phase. A thermodynamically inhomogeneous solution is physically inacceptable 
and has to be amended. The homogeneous free energy should be independent of v. The 
dependence on this parameter is minimized in the saddle point where the free energy f\ does 
not depend on v at least locally. That is, we demand dfx/du = 0. In the homogeneous state 
this equality becomes identity, but when we have an inhomogeneous solution this equation 
just determines an "equilibrium" scaling parameter v eq . 

Comparing free energy (J2J) where parameters q, x and v are determined from the saddle- 
point equations with the Parisi RSB scheme we find equivalence with 1RSB solution 3]. The 



difference between these two constructions lies in their derivation and motivation. While 
Parisi arrives at 1RSB solution from the demand of maximality of the free energy, here we 
derived the saddle-point equations from the effort to restore thermodynamic homogeneity 
(minimize inhomogeneity) . We believe that the latter construction has a more profound 
physical justification. 

The local thermodynamic homogeneity achieved by free energy (J2J) is not enough to 
become a thermodynamically consistent solution. It would be if this newly won solution were 
independent of a new scaling (replicating) of the corresponding spin space. To check whether 
free energy (J2J) is homogeneous, i. e., independent of global scalings of the phase space, we 
write the number of real replicas as a product v = v\Vi, whereby e. g. V\ be determined 
from stationarity equations resulting from Eq. (J2J). We then have two types of real replicas. 
The first one with indices a, b = l,...,z/i determines the 1RSB solution, while the other 
with indices a', b' = 1, . . . , z/ 2 serves to a replication of the phase space in which free energy 
fi was derived. Using again the replica-symmetric ansatz for averaging we obtain only three 
parameters that we denote (Q a ^ b ' a ' a ') aa = q + Xi, {Q ab,a '^ h ') aa = <? + X2, (Q ab,a ' b ') a ^ = q. It 
is again a straightforward task to derive an explicit representation for the free energy in the 
^2-times replicated spin space leading to Eq. (J2J). Not surprisingly we find that the resulting 
free energy formally corresponds to the Parisi 2RSB solution. 

We can proceed with phase-space scalings in each higher hierarchical state with 
q, xi, ■ ■ ■ , Xk and u%, . . . , vk as order parameters. We would end up with the Parisi full 
(infinite) RSB scheme. However, this is not the objective of this construction. We should 
produce at the end a solution that is globally thermodynamically homogeneous, the limit 
n — > is scale independent. 

III. GLOBAL THERMODYNAMIC HOMOGENEITY 

Determining all the physical parameters q, Xi, ■ ■ ■ , Xk and the geometric ones Vi, . . . , 
from the saddle-point equations of the i^RSB free energy, we are left with the number 
of hierarchies K as the only free parameter in this construction. This number will be 
fixed by reaching the global thermodynamic homogeneity. The solution becomes globally 
homogeneous when the free energy does not decay into a higher nontrivial hierarchy, i. e., 
when the number of parameters Xii ■ ■ ■ > Xk does not increase when the phase space (replica 
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FIG. 1: Stability parameters Aj and Ag for a) SK model at T/J = 0.25 and b) 3-spin model at 
T/J = 0.4. Equilibrium v eq is indicated showing that the 1RSB in the SK model in inhomogeneous, 
while in the p-spin model is globally homogeneous. 

index) is scaled. 

Having the 1RSB free energy (J2J) we have two possibilities how a 2RSB solution with 
Xi > X2 > may arise. First, x Xi an d Xi = X2 becomes unstable. It happens if 

Ai = l-/3 2 ((p,(l-* 2 ) 2 ) A ) <0, (3a) 

where we denoted t = tanh \p{r)y/q + A-^/x)] , (X(\))\ = f^VX X(A) and p v = 
cosh^ [P(rjy/q + Xy/x)\ I (cosh^ + Xy/X)])\- Notice that condition (jo^aj) is formally 

equivalent to the smallest replicon for 1RSB. Second, x ~ ¥ Xi an d X2 = becomes unstable, 
which is the case if 

A = 1 - P 2 ( [1 - ( P ut 2 )x + v {( P ut 2 )x - Ml)] 2 ) ri < • (3b) 

Hence free energy (J2J) becomes thermodynamically homogeneous if both numbers Ai and Ao 
from Eqs. Q are nonnegative for the equilibrium scaling parameter u eq . It is never the case 
in the SK model but it may happen above the Gardner temperature in the p-spin infinite- 
range model as illustrated in Fig. In the p-spin model the stability parameters are only 
slightly modified A 1 = 1 - (3 2 p(p - l)ql~ 2 {. . .)„/2, A = 1 - (3 2 p(p - l)q p ~ 2 {. . .) v /2. For 
p = 2 it is the SK, for p > 2 there is always qo = leading to Ao = 1 |4|. 

To conclude, we demonstrated that the demand of thermodynamic homogeneity can be 
build in into the replica trick as scale invariance of the replica index n. Successive scalings of 
the replica index lead to a formal dependence of the solution on the scaling parameter. The 
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replica symmetric solution depends on the scaling parameter virtually only in the spin glass 
phase and hence gets unstable. The RSB scheme of Parisi simultaneously with conditions 
for the global homogeneity (scale invariance) are then derived from successive scalings of 
the replica index by minimizing the incurred inhomogeneity (demanding local homogeneity 
or scale invariance). 
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